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Abstract 

We study the thermalization of a strongly coupled quantum field theory in the presence of 
a chemical potential. More precisely, using the holographic prescription, we calculate non- 
local operators such as two point function, Wilson loop and entanglement entropy in a time- 
dependent background that interpolates between AdS^+i and AdS^+i-Reissner-Nordstrom 

for d = 3, 4. We find that it is the entanglement entropy that thermalizes the latest and 
thus sets a time-scale for equilibration in the field theory. We study the dependence of the 
thermalization time on the probe length and the chemical potential. We find an interesting 
non-monotonic behavior. For a fixed small value of Ti and small values of fi/T the thermal- 
ization time decreases as we increase fi/T, thus the plasma thermalizes faster. For large values 
oi jJb/T the dependence changes and the thermalization time increases with increasing n/T. 
On the other hand, if we increase the value of {Ttj this non-monotonic behavior becomes 
less pronounced and eventually disappears indicating two different regimes for the physics 
of thermalization: non-monotonic dependence of the thermalization time on the chemical 
potential for Tl «1 and monotonic for >> 1. 
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1 Introduction 



A standard method to study near-equilibrium physics is to shghtly perturb the system away from 
the equihbrium state and study the hnear response of the system to these perturbations. One 
of the main virtues of this method is that one can compute observables by calculating correla- 
tion functions in the equilibrium ensemble; The response functions are given in terms of Green's 
function evaluated in the equilibrium state. However, calculating Green's functions in a strongly 
coupled theory is not an easy task. Fortunately, the gauge/gravity duality is a powerful tool to 
compute correlation functions for a large class of strongly coupled theories. It maps the problem to 
a classical gravity calculation in higher dimensional curved spacetime that in the case of conformal 
theories is asymptotically anti-de-Sitter (AdS). The gauge/gravity duality has proven to be very 
useful in the study of near-equilibrium physics (see [I] and [2] and references therein). However, 
out-of-equilibrium processes cannot be understood using linear response theory. The study of out- 
of-equilibrium processes in strongly coupled field theories involves non-trivial temporal dynamics 
and is an extremely challenging problem. The gauge/gravity duality translates this problem into a 
less daunting though equally fascinating one: the study of time-dependent gravitational dynamics, 
and in particular, black hole formation. 

Non-equilibrium states in the gravity dual can be created by turning on time-dependent back- 
ground fields at the boundary. For example if a field is given a temporal dependence of compact 
support (0, 6t) at the boundary it creates a wave that propagates into the bulk. It was shown in 
[3] that, under suitable conditions, the gravitational collapse of this wave will form a black brane 
in the bulk. This process of black hole formation via gravitational collapse is expected to be the 
dual description of the thermalization process in a strongly coupled field theory. 

The goal of the present work is to use the holographic setup explained above to study the 
thermalization in a strongly coupled theory in the presence of chemical potential. A practical mo- 
tivation for this work comes from the fact that the equilibration time observed at the Relativistic- 
Heavy-Ion-CoUider (RHIC) is much shorter than predicted via perturbative approaches to ther- 
malization. This indicates that the strong coupling is crucial in understanding the thermalization 
process of the Quark Gluon Plasma (QGP). Hence, this is a natural arena where gauge/gravity du- 
ality can be useful. Indeed, in |H|5] the authors studied thermalization probes in a time-dependent 
background with vanishing chemical potential. They found that, unlike what is expected preturba- 
tively, the thermalization proceeds top-down i.e. UV modes thermalize first, IR modes thermalize 
later. They also found that among the probes studied it is the entanglement entropy that thermal- 
izes the latest and thus sets the time-scale for thermalization. The thermalization time scales as 
Tcrit ~ ^/2, where i is the typical length of the probe. A numerical estimate yields Tcrit ~ 0.3/m/c. 
Since the RHIC and the Large Hadron Collider (LHC) processes occur at different values of the 
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chemical potential, a natural question to ask is how do these results change when we consider a 
theory with a non-vanishing chemical potential. However tempting, our intention is not to make 
precise numerical predictions based on a toy model but to learn qualitative features related to the 
effect of the chemical potential in the thermalization process of a strongly coupled theory. 

Beyond the possible connection to RHIC and LHC physics, the study of non-local probes 
in a time-dependent background is an interesting problem in its own right and it can teach us 
important lessons about out-of-equilibrium physics. 

To explore the thermalization of a strongly coupled field theory with chemical potential we con- 
sider a time-dependent spacetime that interpolates between AdS^+i and AdS-Reissner-Nordstom 
(AdS-RN) background of the same bulk dimensions. This {d + l)-dim interpolating background, 
also known as the AdS-RN-Vaidya spacetime, can be understood as describing the gravitational 
collapse of a thin shell of charged null dust. We probe the theory with several non-local observ- 
ables: two point functions, Wilson loop and entanglement entropy. We find that, as in the case of 
vanishing chemical potential, the thermalization is top-down and it is the entanglement entropy 
that thermalizes the latest and thus, sets the time scale for thermalization in the field theory. 

Undoubtedly, our most interesting result is that the behavior of the thermalization time as 
a function of the chemical potential seems to separate the physics in two regimes: ^ 1 and 
Ti ^ 1. For a fixed small value of Ti and small values of /i/T the thermalization time decreases 
as we increase /i/T, thus the plasma thermalizes faster. For large values of /i/T the dependence 
changes and the thermalization time increases with increasing /i/T. On the other hand, if we 
increase the value of Ti this non-monotonic behavior becomes less pronounced and eventually 
disappears altogether. Hence, we observe two different regimes for the physics of thermalization: 
non-monotonic dependence of the thermalization time on the chemical potential for Ti <^ 1 and 
monotonic for Ti ^ 1. The characteristics of the non-monotonic behavior depend on the spacetime 
dimension, and is suppressed in lower dimensions. 

The fact that the physics arranges itself to manifest different qualitative features in the above 
two regimes tempts us to label them as follows: one "classical" and the other "quantum" . For a 
system in thermal equilibrium, Ti ^ 1 corresponds to a classical regime and T£ ^ 1 to a quantum 
regime. Somewhere in the middle these two behaviors connect smoothly demonstrating a classical- 
to-quantum transition. We do observe qualitatively different physics as far as the behavior of the 
thermalization time is considered in these two regimes. However, we are dealing with a system 
truly away from equilibrium and it is not clear to us whether a "classical" or a "quantum" regime 
would stand as a precise notion for such processes, and if so, how the intricacies of such regimes 
interplay. 

It should be emphasized that the background we are working with is ab initio "phenomeno- 
^Or, should it be herself! 
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logically motivated" rather than one obtained from rigorously solving Einstein equations to study 
the formation of a black hole in AdS-spacetime. Another approach towards analyzing such ques- 
tions is to start with an out-of-equilibrium initial field configuration and study the thermalization 
process by solving for the gravitational background itself. Valiant efforts have been made along 
such directions in e.g. |6]-[18] by considering colliding gravitational shock waves. It remains to be 
seen whether a generalization of |3] in the presence of a charged matter yields the desired form of 
the AdS-RN-Vaidya background that we work with. 

This paper is organized as follows: In Section 2 we introduce the holographic setup i.e. we 
present the bulk action and the background metric of AdS-Reissner- Nordstrom in {d + 1) dimen- 
sions. Next, in Section 3, we study the equilibrium behavior of the non-local observables mentioned 
above in (i = 3, 4. Sections 4 and 5 contain our main results. In section 4 we present the gener- 
alization of the Vaidya metric to AdS-RN backgrounds by using a charged null dust to form the 
black hole. We then explore the non-equilibrium behavior of the non-local probes and determine 
the thermalization time for each probe. In Section 5 we summarize our results and elaborate 
on possible future directions. A detailed analysis of minimal surfaces and apparent horizons in 
AdS-RN-Vaydia is found in Appendix A. Appendix B gathers some standard facts related to the 
embedding of 4 and 5 dimensional charged AdS black holes in 10 or 11-dim supergravity. 

Note Added: While this paper was close to completion, we became aware of ref. [1^ which 
partially overlaps with our work. 

2 The bulk action and the background 

We will begin with the AdS-RN background and analyze equilibrium properties of the dual CFT by 
looking at two-point correlation function, Wilson loops and entanglement entropy. The holographic 
prescription for computing such observables is to compute minimal surfaces of appropriate dimen- 
sions in AdS-space. For generality we will analyze these observables in an AdS^+i-background, 
where d = 3,4. As far as the background is concerned, our approach is seemingly "phenomeno- 
logical" or the so called "bottom up" where we simply obtain these backgrounds as solutions to 
some effective gravity action in the corresponding dimension. However, for c? = 3 and d = A — 
which are the ones that we will study here — one can embed the corresponding solutions within 
(11 or 10-dimensional type IIB) supergravity. Also note that within supergravity one usually 
obtains a more general multi-charged black hole solution. The backgrounds that we consider here 
corresponds to setting all these charges equalj^ It is also noteworthy that one can obtain such 

^The reduction of type IIB supergravity yields an A/" = 8 gauged supergravity in d = 4 with 50(6) gauge 
group. This admits a consistent M = 2 truncation coupled to two Abelian vector multiplets with an U{lY Cartan 
subgroup of the full 5*0(6). In general this admits a three charge black hole solution. Similarly, the 5 ''-reduction 
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charged black hole solutions in = 6 as well which comes from the S"^- reduction of 1 1-dimensional 
supergravity; we will, however, not consider this case. 

AdS-RN backgrounds are solutions of Einstein-Hilbert action with a -ve cosmological constant 
coupled to a Maxwell field. Let us start with the effective action of the form 

So = Q I d'+'xV^ {R-2K)-\j d'^'x^gF.^F^'^^ , (2.1) 

which gives the following equations of motion 

R^u -\{R- 2A) g,, = g'^PP^.F^, - ^-g,, {F^^ F^^) , (2.2) 



dp [V^g^'g'-'F,,] = . (2.3) 



The solutions to (2.2) and (2.3) for a general d{> 3) are explicitly given below: 



\ ' f{z) J 2L2 ' 

f{z) = l-Mz'+^^^Q'z'^'''K A = Q{zi~'-z'-'). (2.4) 

Here x is a (c? — l)-dimensional vector, M is the mass of the black hole and Q is the charge. The 
constant zh denotes the location of the horizon which is obtained by solving f{z) = 0. In order 
for the one-form At to be well-defined at the horizon, we have arranged At{zH) = 0. The case of 
d = 2 is special and is not captured by the solution given above. In this case, the background 
takes the following form 

ds'^ = — { -f{z)dt^ + ^7-- + dx' 
z^ V f{z) 

f{z) = 1 - Mz^ + ^ log (z/L) , At = Q log (zh/z) , A = --^ . (2.5) 

In the coordinate system used above the boundary is located at 2; = 0. Having a non-zero 
electric field corresponds to having a chemical potential in the dual field theory and we simply 
read off this chemical potential from the vector field At as 

fi = limAtiz) . (2.6) 
Note that the above definition does not work for d = 2|£] 



of 1 1-dimensional supergravity admits an A/" = 2 consistent truncation with U{1)'^ gauge group, which in general 
admits a four charge black hole solution. See e.g. pOj for a discussion of such solutions. Some features of such 
solutions are also discussed in Appendix A. 

■^The case of c? = 2 is rather special. In this case, the identification of the source and the VEV are subtle and 
the chemical potential should be identified as the sub-leading term as 2; — > 0: fi = Q log (zh/L). For more details 
on this issue, see RTl. 
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On the other hand, the temperature of the black hole is given by 



(2.7) 



We will briefly investigate the behavior of the temperature as the chemical potential of the system 
is cranked up. This will be a numerical endeavor and for that purpose we will set L = 1, which 
implies we will be measuring the temperature and the chemical potential in units of L. Typically, 
for a fixed value of M, Q has an upper bound beyond which there is no real positive solution of 
f{z) = As this upper bound is reached, the background temperature approaches zero and we 
approach an extremal solution with zero temperature but finite horizon radius (giving rise to finite 
entropy). The corresponding supergravity solutions (in d = 3,4) are the BPS black holes discussed 
in [22]. These extremal solutions possess a naked singularity. The behavior of the temperature 
with the chemical potential (for fixed M = 1) is shown in fig. [Tj It is noteworthy to comment 
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Figure 1: The dependence of temperature on the chemical potential in (a) d = 3 and (b) d = 4. 
We are measuring the temperature and the chemical potential in units of the AdS-radius. We 
have also set M = 1. 



here that the physics of AdS-RN background has many interesting properties investigated earlier 
in e.g. [231121] and more recently from a more application-towards-condensed matter theory point 
of view reviewed in e.g. |25j . 

Before we proceed further, note that with the following change of coordinate 

dz 

dt = dv + Y ^ (2-8) 
''^This is true for all cases except d — 2 where Q can be as large as possible. 
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the above AdS-RN metrics given in equations (2.4) can be brought in the following fori 



l2 

ds^ = {-fiz)dv'^ - 2dvdz + dx'^) , = AJz) , (2.9) 

^2 J 



where x is a (c? — l)-dimensional vector. The form of the background in (2.9) is suitable for 
expressing the Vaidya metric, which we will discuss later. Note that the form above is generic for 
any d^ where the information of the dimensionality is entirely carried by the function ji^z) and the 
electric field A,,,{z). For our purposes, the explicit expression of A^{z) will not matter at all since 
ultimately we will study minimal surfaces which do not couple to this vector field. 



3 Non-local observables in equilibrium 

Before proceeding to discuss aspects of the thermalization, we begin by exploring the equilibrium 
behavior of the non-local observables in the presence of a chemical potential. The non-local ob- 
servables we will study here are the two-point functions (for operators with large dimensions), 
expectation values for Wilson loops and entanglement entropy. To study this we take the back- 



ground in the form presented in (2.9). 



3.1 Two point function: the geodesic approximation 

The idea here is to study the thermal Wightman function of an operator with large conformal 
dimension. As shown in [26], in this limit the equal-time two point function is given by 

(C»(t,f)C»(t,f' )) ~ e"^^*!^— ' , (3.10) 

where A is the conformal dimension and ^thermal is the renormalized geodesic length. The above 
expression makes use of a saddle point approximation. 

We want to consider space-like geodesic connecting the two boundary points: (t,Xi) = 
(to, — •^/2) and (t',a;'J = (to,-^/2), where (whenever necessary) all other spatial directions are 
identical at the two end points. Such a geodesic is parametrized by f = v{x) and z = z{x) where 
x\ = X. The boundary conditions satisfied by this geodesic is 

z{-l/2) = zo = z{l/2) , v{-l/2) =to = v{i/2) . (3.11) 

Here zq is the IR radial cut-off near the boundary. 
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Note that introducing the Eddington-Finkelstein coordinate in (2.81 yields a gauge field of the form: Ay — A{z) , 



Az = A{z)/f{z), where A(z) is given in (2.4). We can make a gauge choice to recast all physical information by 
choosing a gauge Ay{z) only. 
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The geodesic length is 



j-l/2 J 

C= dx- [1- fv'^ -2v'z'] 
J-e/2 ^ 



1/2 



dx 



(3.12) 



There are two conservation equations: one because the Lagrangian is independent of x and the 
other because the Lagrangian is independent of v. The first one gives 

1 _ 2z'v' - fiz)v'' = 4 , 

where 2* is the value of z[x) at the midpoint. Using the definition of v we get 

dz dv 1 dz 

dx 



V = t, 







f{z) dx 



we can substitute v' in favour of z' in (3.13) and obtain 

v2 



dz 
dx 



1/2 



(3.13) 



(3.14) 



(3.15) 



where the positive sign is taken for x > and the negative sign is taken for a; < 0. The boundary 
separation length is then read off as 

-1/2 

(3.16) 



where Zq is the IR radial cut-offj^ The corresponding integral can be analytically carried out in 
purely AdS-background to yield 



^AdS = 22;* . 



(3.17) 



On the other hand, using the conservation equation, the geodesic length can be obtained to be 

dz 1 



C = 2L 



Z 



(1 - S) 



(3.18) 



Note that both formulae in (3.16) and (3.18) apply for any d. The geodesic length is a divergent 
quantity; we can regularize this length by subtracting off the divergent piece in pure AdS-space. 
In pure AdS-space, this length is given by 



AdS 



2L 



dz 



1 



-2Llog 



-2L\og(j-^ +2Llog 



2£ 



AdS 



(3.19) 



^Thc length integral is perfectly convergent, therefore we can take zq — without any issue. 
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where we have used (3.17) to substitute for in the above expression. We are interested in the 



finite part of this geodesic length. We therefore consider the following renormalized length]^ 



Now we can numerically study how ^thermal behaves in various dimensions as we tune the chemical 
potential. Recall that our underlying theory is conformal, therefore the only relevant parameter 
that we can vary is a dimensionless ratio constructed from the temperature and the chemical 
potential. In {d + 1) bulk space-time dimensions (dual to a theory in d boundary space-time 
dimensions) 

[T] ~ > M ~ v^-^ ■ (3-21) 

^ ^ length length ^ ' 



1 //i 



Thus we can consider 

X,.„ = {^) (3.22) 

to be the relevant parameter that we will vary. In practice we can vary X{d) by first fixing M = 1 
and then just varying the parameter Q till it reaches its maximum value beyond which naked 
singularities appear in the gravitational background. It is straightforward to check that for this 
range of values for Q, X(d) ^ [0, oo\. 

For a given value of X(d)^ we can now study the behavior of Ahermai as a function of the boundary 
separation length. To generate these curves we do the following: To begin with we set L = 1, 
which means the dimensionful quantity Ahcrmai (and later the area or the volume corresponding to 
the Wilson loop or the entanglement entropy calculations) is measured in units of the AdS-radius. 
Now we fix M = 1 and for a given value of Q start with a 2;* very close to the horizon. Next we 
keep changing z^ till we reach — > 0. Each choice of z^ generates an unique value of ^thermal 



and i through equations (3.18) and (3.16) respectively. Finally we plot this result. The behaviour 
of iZthcrmai is showu in fig. |2} The boundary length ^ will be expressed in terms of the boundary 
temperature, T. It is clear from fig. |2| £thermai ~ {Td^) for large T£. 

The general observation we see from these plots is increasing X{d) monotonically increases the 
value of Ahermai for a given boundary separation £ in both = 3, 4. The dependence of i2thcrmai 
on the dimensionless ratio X{d) is shown in fig. [2]^c): the monotonically increasing function is 
generally non-linear — however — for X{d) ^ 1, -Cthcrmai ~ X{d) with a slope that depends on d. 
Non-linearities appear only for small values of X(d)- Also, we have checked explicitly that in the 
linear regime the slope of ^thermal vs X{d) curve depends on the fixed value of {Ttj. This will be a 
generic feature in all the observables we will consider here. 



''Note that in f^, the authors include a finite piece 2 log 2 in the subtracted part. This should not change the 
qualitative behaviour since it's a constant. 
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Figure 2: Left panel (a): The case d = 3. The blue curve corresponds to X(3) ~ 0.003 and the red 
curve corresponds to X(3) ^ 22.4. Right panel (b): The case d = A. The blue curve corresponds 
to X(4) ~ 0.002 and the red curve corresponds to X{4) ~ 18.3. In (c), we have shown how >Cthermai 
scales with X(d) ^or both d = 3 with {A'itT)£ = 0.6 and d = 4 with {AttT)£ = 2. Here ^thermal is 
measured in units of the AdS-radius, L. 



3.2 Space-like Wilson loops 

Wilson loop is another gauge invariant non-local observable that can probe thermal properties of a 
field theory, e.g. the expectation value of the Wilson loop can detect confinement /deconfinement 
transition in theories like QCD. Here we will study the thermalization of space-like Wilson loops 
from a holographic approach. 

In a gauge theory, the Wilson loop operator is defined as a path ordered contour integral over 
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a closed loop C of the gauge field 

W{C) = ^Tr (Ve^c^^ , (3.23) 

where A is the gauge field and is the rank of the gauge group and V denotes path ordering. In 
the AdS/CFT correspondence, the expectation value of the Wilson loop is related to the string 
partition function 

{W{C)) = I I)Se--^(^) , (3.24) 



where S is the string world sheet which extends in the bulk with the boundary condition 9S = C 
and AiTj) corresponds to the Nambu-Goto action for the string. In the strongly coupled limit, we 
can simplify the computation by making a saddle point approximation and evaluating the minimal 
area surface of the classical string with the same boundary condition SSq = C 

{W{C)) = e-'^(^o) , (3.25) 

where Sq represents the minimal area surface. In the AdS/CFT correspondence, such computation 
of Wilson loop operators were first done in [27] . Here we will consider both rectangular and circular 
Wilson loops, which are schematically shown in fig. [3j 

3.2.1 Rectangular Wilson loop 

A rectangular strip Wilson loop can be parametrized by the boundary coordinates {xi,X2} with 
the assumption that this infinite rectangular strip is invariant under the X2-direction. Thus the 
corresponding minimal surface is parametrized by: z{x) and w(x), where x = xi. As in the case 
of geodesies, the boundary conditions are 

z{-i/2) = zo = z{i/2) , v{-i/2) =to = v{i/2) , (3.26) 

where i is the length of the rectangular Wilson loop along the x ^-direction. The Nambu-Goto 
action is given by 

^NG = 7^ — 7 / -^[l-fv -2vz) , 3.27 

where R is the length along the a;2-direction. Since there is no explicit x-dependence in the 
Lagrangian, the corresponding conservation equation is given by 

I- fv'^ -2v'z' = (^)' , (3.28) 
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(a) (b) 

Figure 3: A schematic diagram of the Wilson loops of different shapes and the corresponding 
minimal area surfaces: rectangular in (a) and circular in (b). 



where is the midpoint of z{x). Using the definition of v from (3.14) we can obtain 

dx , 1 r . 4 
dz 



v'2 



/ 



z 



V7L 



Z 



(3.29) 



where the positive sign is taken for a; > and the negative sign is taken for x < 0. The boundary 
separation length is obtained by integrating the above equation from zq to z*. Once again this 
length can be analytically computed for the pure AdS-case to give 



^AdS 



£(|) 



(3.30) 



Using the solution in (3.29) in (3.27), the area of the minimal surface in purely AdS-background 
is given by 

2 



^AdS = Oi'-A 



NG 



TTZo 



fAds Vr(l/4)y' 



(3.31) 
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where the first term in the above expression is divergent. Subtracting this diverging piece from 
the AdS-RN case, we get the following area of the minimal surface 

vAthermal — « A^G — / ^ — , „ , , , , • [6.61) 



Note also that the formula in (3.32) is valid for any d since the information of the dimensionality 
is encoded in the function /. Here also we use a similar technique as outlined before to generate a 
curve of ^thermal (mcasurcd in units of appropriate powers of L) vs Tl. The dependences are shown 
in fig. |4j The general observation is once again for increasing X{d)) -^thermal increases monotonically: 
For large values of T£, we find ^thermal ~ [Ti] with a slope that depends on X(d)- In|4]^c) we have 
shown how the area functional scales with the dimensionless ratio for fixed values of {Tt). 
Here also we find that the generic behavior is non-linear, but it becomes dominantly linear for 
large values of X{d) with a dimension-dependent slope. Non-linearities only show up for small 
values of X{d)- 

3.2.2 Circular Wilson loop 

Let us now discuss the circular Wilson loop. At the boundary we choose a 2-dimensional plane 
{xi,X2\ and rewrite it in the polar coordinate {p, 0} 

dx\ + dx\ = dp^ + p^dcjP' . (3.33) 

Using the azimuthal symmetry in the 0-direction, the minimal area surface can be represented by 
z{p) and v{p). The Nambu-Goto action for this Wilson loop is given by 

a Jo ^ 

where ' = d/dp and we have integrated over the azimuthal angular direction that yields a factor 
of (27r). Using the definition of v once again we can write this as 



Note that in this case we do not have any conservation equation. The equation of motion for z{p) 



obtained from (3.35) is given by 



,o 1 /2 1 df\ z' 2f , , 

z" + z'^— + z'U--—-^)+- + ^ = 0. 3.36 
pf \z 2fdzJ p z 

The boundary conditions we should use are given by 

z{0) = z, , z'{0) = . (3.37) 
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Figure 4: Left panel (a): The case d = 3. The blue curve corresponds to X(3) ^ 0.003 and the red 
curve corresponds to X(3) ^ 22.5. Right panel (b): The case d = A. The blue curve corresponds 
to X(4) ^ 0.002 and the red curve corresponds to X(4) ~ 18.3. In (c), we have shown how ^thermal 
scales with X{d) for both d = 3 and d = 4 for fixed [^■kT)1 = 2. ^thermal has been measured in 
units of the AdS-radius. 



By expanding near the p = point, we get 



/I - Mzl Q% » 2 / N 

^middle(p) = Z^- \ h —Z^ ] p + ... . (3.38) 



2z, 4 

We use the above expansion at p = e to impose the boundary conditions and solve the equation 



(3.36) using Mathematica's NDSolve, where e is some small numberjj To generate the desired 
curve we now have to carry out the following steps: First we fix M (set to unity) and Q. Then 



^Typically we have used e ^ 10 
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we choose some z^, which is numerically close to the event-horizon. For these choices of the input 
parameters and the above boundary conditions we first solve the differential equation numerically. 
To read off i?, we now impose the (radial) IR boundary condition: z{R) = zq, where Zq is the IR 
cut-off. Using this numerical solution we can generate an unique value for v4thermai for a given R. 
Now we can vary and keep repeating the same process to generate the curve we want. 

Before presenting the results, let us recall again that this area is formally a divergent quantity. 
To extract the divergent piece, we can again just focus on the pure AdS-case. In this case (which 
corresponds to setting M = and Q = 0), a. simple solution of (3.36) is given bj|^ 



z{p) = ^/W^ , (3.39) 

which ultimately gives 

R 

AAdS = — 7 , 3.40 

a' Zq a' 

where zq is the radial IR cut-off. Thus the finite contribution to the thermal Wilson loop is 
obtained to be 

Ahermal = « ^NG • (3-41) 

Zq 

The numerical results are shown in fig. [5j It is clear that the expectation value of the Wilson loop 
depends on the dimensionless parameter X(d)- Our general observation again seems to hold here, 
i.e. increasing X(d) increases the area of the minimal surface for a given value of [TR). Note that 
in this case, for large values of (TR), ^thermal ~ (TR)"^. 



3.3 Entanglement Entropy 

Consider a quantum field theory with many degrees of freedom and assume that the zero tem- 
perature ground state of the system is described by the pure ground state denoted by |\I'), 
which does not have any degeneracy. The von Neuman entropy of this system, defined as 
•S'totai = -trplogp = 0, where p= |^) (^|. 

Now we can imagine dividing the system into two subsystems A and B. The total Hilbert 
space now factorizes as "Htotai = 'Ha ^'Hb- Now let us imagine an observer who has access only to 
the subsystem A, consequently the relevant density matrix for this observer is the reduced density 
matrix defined as 

PA = trBp. (3.42) 

^This solution was obtained in [28 by taking an infinite straight Wilson line and then applying conformal 
transformation to map the line to a circle. 



14 



■^thermal 

1.0 
0.5 

0.0^ 
-0.5 

-1.0 



^ (AnT) R 



(a) 



'^thermal 

0.5^ 



0.0 



-0.5 



■1.0 



4 6 / 



(b) 



/ 



/ n / (4;rT) R 



Figure 5: Left panel: The case (i = 3. The blue curve corresponds to Xiz) ~ 0.003 and the red 
curve corresponds to X(3) ~ 0.19. Right panel: The case d = 4. The blue curve corresponds to 
X(4) ~ 0.002 and the red curve corresponds to X(4) ~ 0.24. ^thermal is measured in units of the 
AdS-radius. 



The entanglement entropy is now defined as the von Neuman entropy using this reduced density 
matrix 

Sa = -tTAPA^OgpA ■ (3.43) 

Entanglement entropy measures the quantum entanglement between the two subsystems A and B 
and it is non-zero even at zero temperature. Entanglement entropy is also an useful order parame- 
ter in quantum phase transitions. Here we will eventually use this observable to probe the physics 
of thermalization in the presence of a chemical potential and we will explicitly demonstrate that 
this is the observable that sets the time-scale of equilibration since of all the non-local operators 
considered, it thermalizes the latest. 

A prescription for computing entanglement entropy using the AdS/CFT correspondence was 
suggested by [29] and has been analyzed in details in the literature since then. For time-dependent 
background, the covariant proposal for entanglement entropy was suggested in [30]. According 
to this proposal one needs to consider extremal surfaces instead of minimal area surfaces For a 
recent review, see e.g. [STJ- Let us consider bulk AdS^+i space-time. Then the formula for the 
entanglement entropy is given by 

_ Area(7A) , . 

where G^^~^^^ is the {d + l)-dimensional Newton's constant; 7^ denotes the {d — l)-dimensional 
minimal surface whose boundary coincide with the boundary of the region A: d'jA = dA. 

In d = 3, the computation of the entanglement entropy is identical to the Wilson loop compu- 
tation. Thus here we will only consider the case d = 4. As before, we will consider two geometric 
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shapes: infinite rectangular strip or the straight belt and the spherical region. Fig. [3] again serves 
as a pictorial representation of the particular shape considered and the corresponding minimal sur- 
face. Thermalization of the entanglement entropy in d = 2 was first studied in [32] for vanishing 
chemical potential. 



3.3.1 The straight belt 



The straight belt can be parametrized by the boundary coordinates {xi,X2,X3} with the assump- 
tion that this infinite strip is invariant under both X2 and Xs-directions. The corresponding minimal 
surface can be represented by z{x) and v{x) where x = Xi. As before, the boundary conditions 
imposed are 



z{i/2) = zo = z{-i/2) , v{i/2) =to = v{-i/2) . 



The volume functional is given by 

L^dx 



V = A 



(1 - fv'^ - 2v'z') 



1/2 



where V = Area(7yi) . 



(3.45) 



(3.46) 



-1/2 ^ 

Here A is the area that results from integrating over the x"^ and the x^-directions. The equation 
of motion resulting from the conservation equation is given by 

z 



fv'^ - 2v'z' 



(3.47) 



Using the definition of v from (3.14) we can obtain 



1 + 



/ 



dz 
dx 



1 



1/2 



(3.48) 



where the positive sign is taken for a; > and the negative sign is taken for a; < 0. We can 
analytically obtain the length for pure AdS-background to give 

T(2/3) 



r(i/6) • 



Finally the volume functional is given by 



V = 2AL^ 



dz 1 



z'y/fy/l-{z/z,y 



(3.49) 



(3.50) 



which is again a formally divergent quantity. In pure AdS, this volume functional is obtained to 
be 



V 



AdS 



2AL^ 



dz 



1 



AL' 



20 



+ AL' 



;V^ r(-l/3) 

3^2 r(i/6) 



AL^ AL^7i^r{-l/3){r{2/3)y 



+ 



-AdS 



(3.51) 
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Thus the finite part of the volume functional can be obtained to be 

dz 1 1 AL^ 



thermal 



(3.52) 



The behaviour of the entanglement entropy is demonstrated in fig. |6| Here we also notice that 
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Figure 6: The case d = 4. Left panel: The blue curve corresponds to X(4) ~ 0.002 and the 
red curve corresponds to X(4) ~ 18.3. Right panel: the scaling of Vthermai with X{4) for fixed 
(47rr)^ = 1.5. We measure Vthcrmai in units of the AdS-radius and also set A = 1. 



for large values of (T£), Vthcrmai ~ (Tt) with a slope that depends on X{i)- On the other hand, the 
scaling of Vthcrmai with X(4) displays the non-linear and monotonically increasing behavior where 
the non-linearities are washed out for large values of X(4)- Although we have not displayed it in 
the figure, the slope of the linear behavior of Vthcrmai with X{a) for large X(4) depends on the value 
of (T£). 



3.3.2 The spherical region 

Now we consider the case where one of the subsystem is a circular disc. To parametrize this 
circular disc in polar coordinate, we rewrite 

^ dxl = dp^ + p^dnl_^ . (3.53) 
1=1 
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The minimal area surface can now be parametrized by z{p) and v{p). The volume element in this 
case is given by 



'0 

4vr /'dp^fl + ^l . (3.54) 



V = An [ dp^ (1 - fv'^ - 2z'v') 
Jo z 

We do not have the integral of motion anymore and the full equations of motions are obtained 
by varying the functional in (3.54). To avoid clutter, we do not present the specific form of the 
equation of motion here. The boundary conditions are once again 

z{0) = z^ , z\0) = . (3.55) 

As before, solving the equation near the p = region we get 

3-3M4 + 2g^4 2 ^ 
z{p) = z^ p + ... . (3.56) 

bz^ 

In practice we use the above expansion at p = e to impose the boundary conditions. Here e is a 
small number typically of the order of 10^^. 

To determine the divergent piece in the volume, let us evaluate it in the pure AdS-case. The 
solution of the minimal surface is simple: z"^ = — p^ . So we get 

rp{zo) 2 /d2 z \ 

VAdS = 47ri?L^ / dp 2 = 2?^^^ ^ + log — ^ + finite . (3.57) 



Hence the finite part we are interested in is given by 

« p2 / ^,2x1/2 /^2 







z-" 



Vthermal = AtiL' / rfp^ ( 1 + y J - 27rL^ + log . (3.58) 



The dependence is shown in fig. [7| Once again we observe the behavior that increasing X(i) 
increases the volume functional for a spherical subsystem of a given radius. For large values of 
{TR), we recover a quadratic behavior of the volume functional, which is expected from the general 
area-law scaling of entanglement entropy. 



4 Non-equilibrium physics 

4.1 The bulk action and the background 

Our goal here is to study the same set of non-local observables in a time-dependent background. 
This time-dependent background should capture the physics of the formation of a black hole — 
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Figure 7: The case d = A. The blue curve corresponds to X(4) ~ 0.002 and the red curve 
corresponds to X(4) ~ 0.24. The dimensionful quantity Vthermai is measured in units of the 
AdS-radius. 



in the present context — a black hole with a definite mass and charge. To that end, we will use 
a generalized version of the AdS-Vaidya background including a charge for the black hole. For 
obvious reasons, we will call this the AdS-RN-Vaidya background. 



To find the corresponding Vaidya background, we have to couple the above action in (2.1 ) with 
an external source 

S = So + nS^^t , (4.59) 

where k is a constant and we do not specify the form of S^xt- The equations of motion in this case 
will take the following form 

R^.u 'liR- 2A) 9^.. - g'^'Fp^F^, + ^-g^, [F'^^ F^p) = 2 (sTrCSJ+'^ft:) T; 



next 



(4.60) 
(4.61) 



Let us start with the case where the matter field (corresponding to S'ext) is neutral and subsequently 
the black hole formed will be the AdS-Schwarzschild. The metric is the ((i+l)-dimensional infalling 
shell geometry described in the Poincare patch by 

L2 



{^—f{v,z)dv^ — 2dvdz + dx^) , f{z,v) = l — m{v)z" 



(4.62) 



and is known as the AdS- Vaidya background. Here m{v) is a function that captures the information 
of the black hole formation. On physical ground, m(f ) should interpolate between zero (in the limit 
V — )■ — oo corresponding to pure AdS) and a constant value (in the limit w — oo corresponding to 
AdS-Sch). A choice of such a function is 

. . M 



m(v} 



, 1 + tanh — 
2 V vo 



(4.63) 
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Here vq is a parameter that denotes the thickness of the shell. With this choice, the external 
source must yield the following energy-momentum tensor 



2 SttG 



N 



~ 2 dv 



(4.64) 



If we identify = 6^^, then we get [32 



next 



with A;^ = , 



(4.65) 



which is characteristic of null dust. Thus the formation of the black hole is realized by a shell of 
infalling null dust. In [H |5], this metric has been used to study aspects of thermalization. 

Our first goal here is to generalize the Vaidya metric for the AdS-RN background in {d + 1)- 
dimensions. It is clear that now we need a charged null dust for the formation of the black hole. 
As before, we will present the corresponding Vaidya metric for specific values of d^ 

In (i = 2 we get 

L2 



ds^ = (—f(z,v)dv'^ — 2dvdz + dx'^) , = qiv) \ogz 

q{v) 



f{z,v) = l-m{v)z^ + ^j^z^logz , A = - — 



2«:T;f 



dm 2 dq 

^-^log(.)g(.)^ 



'Jfj.v'-'uv ; ""-'ext 



1 dq 
L dv 



5^' 



(4.66) 
(4.67) 
(4.68) 



In (i = 3 we get 



ds^ 



[—f{z,v)dv'^ — 2dvdz + dx'^) 



q{y)z 



f{z, v) = 1 — m{v)z^ + 



q{vf_ 
2L2 



A 



^■^ext — 7 

dv 



2/cT-* = 



dm 
dv 



3 

z , ,dq 
dv 



(4.69) 
(4.70) 
(4.71) 



And, finally in d = 4 we get 

L2 



ds' 



f{z,v) 



(^—f(z, v)dv'^ — 2dvdz + dx'^) , = q{v)z'^ 



m{v)z'^ + 



2q{vf 



^Jext — 21/ — (5^^ 

dv 



2kT: 



cxt 



Z*" , A 

3 o dm 2z 

-Z -r- 

2 dv 



6 

Z2 ' 



dq 



L2 



(4.72) 
(4.73) 
(4.74) 



""^"Note that here we are not careful about the regularity of the one-form at any particular point since the 
horizon is created only in the v ^ oo limit. Perhaps a better way to write the background is to write the field 
strength instead of the gauge field itself. For all our purposes though, this gauge field does not play any role. 
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As before, we will only consider the cases d = 3, 4j^ For all these metrics, we can choose the profile 



for m{v) as given in (4.63) and we are also free to pick an analogous profile for q{v). For future 
purposes, we will set the radius of the AdS-space, L = 1, and thus all dimensionful quantities (the 
length, area or the volume functional as the case may be) will be measured in units of this scale. 



Now we will work with the backgrounds given in (4.69) and (4.72). As previously stated we 
will work with the following mass function 

m{v) = — 1 + tanh — . (4.75) 



2 V ^^0 

To pick a charge function, we note the following [33]: For d = 3 we can rewrite the function f{z, v) 
as follows 

f-i-(^r^f(^)\ (4.76) 



zh{v) J 2 \zh(v 



which gives 

1 



^ m(v) ' 



q{vf = Q''m{vf'^ . (4.77) 



For = 4 we can rewrite the function f{z^ v) as follows 

^zh{v)J 3 \zh{v) 

which gives 



/ = l-(^) +?(^) > (4-78) 



miv) 



q{vf = Q^m{vf'^ . (4.79) 



From now on, we also set M = 1. To appeal to the visual cortex, let us plot how the interpolating 
function m{v) and q{v) behave in various dimensions in fig. [sjand fig. |9} It is clear that in the 
thin shell limit, characterized by fo — ?■ 0, the variation of the mass and the charge functions are 
sharply changing around f = and approximates a step function behavior. In the remainder of 
this paper, we will analyze the observables in these dynamic (time-dependent) backgrounds in the 
thin-shell approximation. 

Before proceeding to the analysis of the probes of thermalization, a notational comment is 
in order: In the previous sections, we have denoted the regularized length, area or the volume 
functional with a subscript "thermal". From here on, we will keep the subscript "thermal" in 
direct comparison with the equilibrium cases but the quantities are not be interpreted as thermal 
ones, rather they have explicit time-evolution. 



A general form of this background in (d + l)-dimensions in given in equation ( A.l ). 
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(a) (b) 

Figure 8: The interpolating functions m(v) and q{v) in d = 3 for various values of vq 
0.01,0.1,0.2 (red, pink and blue). We have set Q = 0.5. 
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Figure 9: The interpolating functions m{v) and q{v) in d 
0.01,0.1,0.2 (red, pink and blue). We have set Q = 0.5. 



4 for various values of vq 



4.2 Spacelike geodesies 

We start by analyzing the two-point function. We consider geodesies with a boundary separation 
along Xi = x-direction (all other spatial directions at both the end-points are the same). The 
profile is described by two functions z{x) and v{x). The length element is 




(4.80) 
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The conservation equation is given by 

1- fv'^ -2v'z' = (^^y . (4.81) 

Using this conservation equation, the two equations of motion corresponding to the variation of 
z{x) and v{x) are obtained to be 

zi," + 2z'v' - 1 + v" (f - ^z^'] = , (4.82) 

.'' + K + |.V+l|.- = o. («3) 

We solve the above two differential equations subject to the following boundary conditions 

z{e) = z^ , z'{e) = + corrections , f (e) = , f '(e) = + corrections , (4.84) 

where e is a small number. In practice, we fix the slopes z' and v' at x = e from the equations 
of motion themselves, which gives the "corrections" terms above. So far and v^, are two free 
parameters that generate the numerical solutions for z{x) and v{x). The boundary data can be 
obtained from this numerical solution 

z{±i/2) = zo , v{±i/2) = t . (4.85) 

Here Zq is the radial IR cut-off and t is the boundary time. To generate the corresponding 
thermalization curves, we do the following: for a fixed value of z^, we keep varying v^, till the 
read-off value of zq is sufficiently low. This generates one desired profile for z{x) and f (x), which 
we can use to compute the length functional. Now we vary z^, and repeat the process again. 



Using the data obtained as explained above we show how the thermalization occurs in fig. 10 



It should be emphasized at this point that we are not being careful about the units in which 



we measure the boundary time in this figure. In what we have shown in fig. [10} the boundary 
time is measured in units of the black hole mass, M, which strictly speaking does not have a 
simple interpretation in terms of any physical quantity of the boundary theory. Moreover, the 
different curves corresponding to different values of X{d) are not obtained for the same equilibrium 
temperature. Thus fig. [TO] and all such subsequent ones in the later sections should be viewed as 
schematic representation of the physics that is going on. 

We need to define a "thermalization time" in order to investigate how it depends on the length 
scale for various values of X{d)- Following [1], we define two different time scales denoted by ri/2 
and Tcrit respectively: 

(i) ri/2 is defined as the time required for the curves to reach half of their equilibrium value. 
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Figure 10: Left panel: The case d = 3. The blue and the red curve corresponds to X(3) ~ 
0.003, 4.45 respectively. Right panel: The case (i = 4. The blue and the red curve corresponds to 
X(4) ~ 0.002, 0.4 respectively. >Cthermai is measured in units of the AdS-radius and the boundary 
times measured in units of the black hole mass M. All curves here correspond to different values 
of (T£), which we do not specify here. 



(ii) Tcrit is defined as the critical time at which the geodesic grazes the middle of the shell at f = 0. 
This is simply given by the following formula: 

Tent = 7^ , 4.86 

where, as before, zq is the IR radial cut-off and determines the value of the boundary separation. 
Once again, we need to specify the scale in which we measure these thermalization times. From 
the point of view of the boundary theory the only relevant dimensionful quantity is the equilibrium 
temperature. Thus we will always study the behavior of the dimensionless quantities (Tri/2) or 
(Trent). 

Now we can investigate how these thermalization times depend on the length-scale i for various 



fixed values of X{d)- The results for T1/2 are shown in fig. 11 For small i, T1/2 grows linearly with 
the property that T1/2 < i/2. The deviation from linearity for large values of i is clear. For 
increasing X{d) this deviation occurs either in the opposite direction or occurs more slowly and 
for large value of i increasing X(d) increases the value of T1/2. For small values of i, T1/2 can 
decrease with increasing X{d)- Thus we already observe possibly interesting and different physics 
dominating two different regimes: (Ti) <^ 1 and (Ti) ^ 1 for increasing values of X{d)- It should 
be noted however that for very small values of i, T1/2 is not very sensitive to the value of X{d)- 

On the other hand, the dependence of Tent with the boundary separation length is shown in 
fig. 12 The qualitative behavior of this thermalization time is much like the one observed for ti/2'- 
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Figure 11: Left panel: The case d = 3. The blue and the red curve corresponds to X(3) ^ 
0.003, 4.45 respectively. Right panel: The case d = 4. The blue and the red curve corresponds 
to X(4) ~ 0.001, 0.3 respectively. 



For small boundary separation, the thermalization time is not very sensitive to the parameter X{d) 
— although there is a range within which increasing X{d) actually decreases Tcrit- However, for 
large values of Tent clearly increases with increasing X(d)- The deviation from linearity for large 
values of i is also quite interesting. For small values of X{d) this deviation from linearity results in 
Tcrit < ^/2; whereas for larger values of X{d) this deviation results in Tent > ^/2, which is seen from 



fig. [12] as the blue and the red curves bending away from each other as i increases. This bending 
away effect is more pronounced for d = 3 as compared to d = 4. 

We can also investigate the dependence of Tent with X{d)- To this end, let us define 

Tcrit= hm Tcrit (4.87) 



and normalize the measured thermalization time in units of r^%. In fig. 13 we have shown the 



results. Once again the results display an interesting interplay of physics for small and large values 



of X{d)- This effect is relatively milder in c? = 3 compared to c? = 4. From fig. 13|^b) we observe 
that for fixed values of (T£), different physics dominates two different regimes: small X(4) and 
large X(4)- Initially the curve decreases before turning back and increasing monotonically. This 
monotonically incasing behavior seems to be a linear one, with the slope depending on the fixed 
value of (T£) and also the dimension we are in. This growth seems to be unbounded, implying 
that for infinitely strong chemical potential, it takes infinitely longer to reach thermalization. On 
the other hand. Tent — measured in units of r^j,;^ — has a minima and both the location and the 
minimum value depends on the fixed value of (T£). From what we have shown in fig. [I3|(b), we 
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Figure 12: Left panel (a): The case d = 3. The blue and the red curve corresponds to X(3) ^ 
0.003, 4.45 respectively. Right panel (b): The case d = 4. The blue and the red curve corresponds 
to X(4) ~ 0.002, 1.1 respectively. 
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Figure 13: Left panel (a): The case of d = 3. Right panel (b): The case of d = 4. Here T^^■^^^ 
denotes the value of Tcrit evaluated at X{d) = 0- 



have about 4%-reduction in thermalization for (Ti) = 0.7/(47r). We will observe that all these 
features are very generic for all non-local probes that we consider in the subsequent sections. 

4.3 Wilson loops 

Once again we will consider two geometric shapes: the rectangular and the circular one and we 
appeal to fig. |3] for a schematic representation. 
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4.3.1 Rectangular strip 

As before the area functional for the rectangular strip is given by 

R /"^/^ dx „ ,o ^ , ,x 1/2 



A 



27r ./_^/2 



{1- fv'^ -2v'z'y 



The conservation equation gives 



12 



2v'z' 



(4i 



(4.89) 



Using the conservation equation and extremizing the area functional we get the following two 
equations of motion 



II . II _c , I I" J , ^ 12" J n 

1 df 



df 



zv" + 4z'v' -2 + v" [2f 



2 dz 



(4.90) 
(4.91) 



Note that compared to the computation of the geodesies, only the second equation changes. We 



again use the same boundary conditions as outlined in (4.84). Some of the resulting plots are 



shown in fig. 14 The behavior of the thermalization time with the length of the Wilson loop 
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Figure 14: Left panel: The case d = 3. The blue and curves correspond to X(3) ~ 0.003, 4.47 
respectively. Right panel: The case d = 4. The blue and red curves correspond to X(3) ^ 
0.002, 0.4 respectively. We measure .4thermai ii^ units of the AdS-radius and the boundary time 
t in units of the black hole mass. 



operator has been displayed in fig. 15 The general behavior of ti/2 is sub-linear. The deviation 



from linearity of ri/2 decreases with increasing X{d) and for large enough i, the thermalization time 
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Figure 15: Left panel: The case d = 3. The blue and red curves correspond to X{3) ^ 0.003, 4.47 
respectively. Right panel: The case d = 4. The blue and the red curve corresponds to X(4) ^ 
0.002, 0.4 respectively. 



increases with increasing value of X{d)- From fig. 15 ^b), we can clearly identify a regime of (Ti) 
where increasing X{4.) actually decreases ti/2- 

On the other hand, the behavior of Tent with the length of the Wilson loop operator is shown 



in fig. 16 It turns out that Tcru exceeds the linear relation in i/2 and for large enough i, the red 



and the blue curves bend away from each other. So increasing X{d) enhances the departure from 



the linear behavior for Tcrit- From fig. 16 b) we can also identify a regime where Tent decreases 



with increasing chemical potential and hence the thermalization happens faster. 



In fig. 17, we demonstrate how the thermalization time depends on X(d) for both d = 3 and d = 
4. As in the case of 2-pt function, this response consists of two different regimes: for small values 
of X(d), thermalization seems faster and for larger values of X{d) thermalization becomes slower 
and approaches infinitely large values for infinitely strong chemical potential. The percentage 
reduction in thermalization time for small chemical potential is about 4% for d = 4 where this 
effect is more pronounced. 

4.3.2 Circular Wilson loop 

Let us now consider circular Wilson loops in the AdS-RN-Vaidya background. The corresponding 
minimal surface is parametrized as before and the area functional has the following form 

^ P 



A{t, R) = / dp^^/l - fv'^ - 2z'v' , (4.92) 
^ 
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(a) (b) 

Figure 16: Left panel: The case d = 3. The blue and red curves correspond to X{3) ~ 0.003, 4.47 
respectively. Right panel: The case d = 4. The blue and the red curve corresponds to X(4) ~ 
0.002, 1.1 respectively. 
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Figure 17: Left panel (a): The case d = 3 for (47rT)^ = 0.4. Right panel (b): The case d = 4 for 
{AttT)£ = 0.7. 



where ' = d/dp. Here the mass and the charge functions are time-dependent and are given by 
(4.75) and (4.77) or (4.79) for d = 3 or = 4 respectively. As before, the equations motion 



resulting from the area functional are messy and therefore we do not present them explicitly. 

To impose the boundary conditions, we can use the rotational symmetry along 0-direction and 
impose 



z(e) = + corrections 



V e 



+ corrections 



(4.93) 
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where and are the two free parameters, which we will describe how to fix. The "corrections" 
in the above is obtained in the following manner: we expand the equations motion near p = e, 
where e is a small number typically of the order of 10^^. Using these expansions we can fix the 
"corrections" as we did in, e.g. (3.37) and ( 3.38| ) f^ Note that this expansion also allows us to set 
z'{e) and f'(e). 

Once these boundary conditions are fixed, we can obtain a numerical solution for z{p) and v{p) 
for every value of z^ and . From this solution we can read off the boundary data 



z{R) = zo , v{R) = t 



(4.94) 



where R is the radius of the circular Wilson loop, Zq is the radial IR cut-off and t is the boundary 
time. In practice we do the following: for a given i?, we fix z^, and keep varying v^, till the obtained 
value of the radial IR cut-off, denoted by zq, is small enough. Then we vary z^ and repeat the 
process. Ultimately this generates the required data to produce the thermalization curves for the 



Wilson loop operator. Some of these curves have been shown in fig. 18 
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Figure 18: Left panel: The case c? = 3. The blue and the red curves correspond to X(3) ~ 
0.003, 0.33 respectively. Right panel: The case d = 4. The blue and red curves correspond to 
X(3) ~ 0.002, 1.1 respectively. We measure ^thermal in units of the AdS-radius and the boundary 
time in units of the black hole mass. The diameter, also measured in units of the black hole 
mass, is fixed to D = 4. 



To extract the behavior of the thermalization time ri/2, we repeat this process for various 
values of the radius and finally obtain the fig. 19 For small values of the behavior of ri/2 is 
primarily sub-linear. The linear regime for small D has a smaller slope in d = 4 as compared to 

^^The only difference is we now have to set the boundary conditions for two functions z(^p) and w(p). 
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the one in d = 3. Increasing X{d) suppresses the sub-hnear behavior of ri/2 and for large enough D, 



ri/2 increases for increasing X{d)- For small values of D, from fig. 19 b) we can definitely identify 
the regime where increasing chemical potential leads to faster thermalization. Also, it is clear 
from the d = 3 case that for moderate values of X(3)5 ti/2 does not change much unless we go to 
very high values of D. 
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Figure 19: Left panel: The case d = 3. The blue and red curves correspond to X(3) ~ 0.003, 0.33 
respectively. Right panel: The case d = 4. The blue and the red curve corresponds to X(4) ^ 
0.002, 1.1 respectively. 



On the other hand the thermalization time, denoted by Tcrit, is shown in fig. |20j For small 
values of X{d), for a relatively large regime of values for D, Tent behaves linearly with a slope very 
close to 1/2. Increasing Xid), promotes a deviation from this linearity for larger values of D and 
the corresponding curve bends away. This means for large enough D, increasing Xid) increases 

^crit • 

4.4 Entanglement entropy 
4.4.1 The rectangular belt 

Let us first begin with the rectangular strip geometry. The volume functional corresponding 
to the minimal area surface can be parametrized by z{x) and v{x), where x = xi, with the 
assumption that this minimal area surface is invariant under the other two planar directions. The 
corresponding volume functional is now given by 

/•^/^ dx 1 /9 

V= ^(i-fy'^-2v'z'Y^\ (4.95) 
J-e/2 z 
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Figure 20: Left panel: The case d = 3. The blue and red curves correspond to X{3) ~ 0.003, 4.47 
respectively. Right panel: The case d = 4. The blue and the red curve corresponds to X(4) ^ 
0.002, 1.1 respectively. 



where ' = d/dx. This gives the conservation equation 



1 £ 12 nil f ^* 

I — jv —2vz=[ — 
V z 



(4.96) 



The two equations of motion resulting from the variation of the volume functional can be written 
as 



zv" + 6v'z' -3 + v'^3f 



1 dl 
2^dv 







(4.97) 
(4.98) 



Once again only the second equation above changes. We will again solve these two equations 



subject to the boundary conditions outlined in (4.84) and using similar method as outhned earlier. 



Some of the curves obtained in this way are shown in fig. 21 



In fig. 22 we have shown how the two thermalization times ti/2 and Tcrit depend on the length 



i of the rectangular region. For small values of X(4), ti/2 behaves linearly with i/2 for small values 
of i: the corresponding slope of this linear behavior is unity. For larger length, ti/2 becomes 
sub-linear. On the other hand, for larger values of X(4) the deviation from linear behavior of ti/2 
is suppressed and for larger length, increasing X{i) increases ri/2- 

The qualitative behavior of T^nt is similar to that of ri/2- However, the slope of the linear 
regime is quite different. For Tcrit, in the linear regime, Tcrit ~ ^ and therefore has almost twice the 
slope compared to the curves for ri/2- Once again we observe that for large enough £, larger X(4:) 
leads to larger Tcrit- Finally, in[22|^c) we have shown the dependence of Tcrit with the dimensionless 
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Figure 21: The case of d = 4. The blue and the red curve corresponds to X(4) ~ 0.002, 1.1 
respectively. Vthermai is measured in units of the AdS-radius. Here we have fixed £ = 2 in units of 
the black hole mass. We do see the appearance of the swallow-tail behavior in the thermalization 
curve for entanglement entropy. 



ratio X{4); which gives a similar result as we have observed earlier, specifically we do seem to have 
two different regimes for X{4) distinguished by a faster or a slower thermalization corresponding 
to small and large values of X(d)- A similar behavior is observed as (T£) is varied for fixed values 

of X(4). 



4.4.2 Spherical region 

Now we will investigate the case of a spherical region. To this end, we parametrize the spherical 



disc in polar coordinates as in (3.53) and represent the corresponding minimal are surface by z(p) 



and v{p). The corresponding volume functional is given by 



R 



P 



V = Att dp^-^l - fv'^ - 2v'z' , 



(4.99) 



where ' = d/dp. We follow similar methods as outlined in section 4.3.2 in obtaining the thermal- 



ization curves and a few representative curves are shown in fig. 23 



We have demonstrated the behavior of the thermalization time with the diameter of the spher- 



ical region in fig. 24 In this case, ti/2 has a sub-linear behavior for the range of the diameter we 
have explored for both small and relatively large values of X(4)- The linear regime grows slower 
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Figure 22: The case of (i = 4. Left panel (a): rx/2 ^-s a function of the boundary separation i. 
The blue and the red curve corresponds to X(4) ^ 0.002, 1.1 respectively. Right panel (b): Tcrit 
as a function of the boundary separation I. The blue and the red curve corresponds to X(4) ^ 
0.002, 1.1 respectively. In (c) we have shown liow T^rit behaves with X(4) for fixed (47rr)£ = 0.5. 



than D jl. Also, for the range of the diameter that we have explored, increasing X(4) decreases 

n/2- 

The behavior of Tent, on the other hand, is in agreement with previously obtained qualitative 
features. For small value of X(4); ^crit ^ D/2 and this linear behavior persists till larger values of 
D. However, as X{4) is increased, Tcrit exceeds this linear behavior and the corresponding curve 
bends away. 
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Figure 23: The case of d = 4. The blue and the red curve corresponds to X{4) ~ 0.002, 1.1 
respectively for fixed diameter Z) = 4 (in units of the black hole mass). Vthermai is measured in 
units of the AdS-radius. 
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Figure 24: The case of d = 4. The blue and the red curve corresponds to X(4) ^ 0.002, 1.1 
respectively. 



5 Summary and Discussion 

We have explored, in details, the behavior of thermalization time as the chemical potential is 
varied by probing the following various non-local observables: two-point function, Wilson loop 



35 



and entanglement entropy. In this section, we briefly summarize and review some crucial features 
of the observations we made. For simplicity, we will only comment on the behavior of r^rit, but 
the behavior of ri/2 is qualitatively similar. 

At vanishing chemical potential, which in our notation would be denoted by X{d) = 0, we 
reproduce the results obtained in |ll|5]. The generic features in the absence of a chemical potential 
are the following: for the two-point function, the thermalization time Tent = i/2 in AdSa but Tent < 
£/2 in AdS4^5 and deviates from linearity. The Tent = i/2 behavior is understood from the dual 
(l+l)-dim CFT. On the other hand, the deviation from linearity in higher dimensions is interpreted 
as an indicator of a faster-than-causal thermalization possibly resulting from the homogeneity of 
the initial configuration In higher dimensions, the Wilson loop or the entanglement entropy 
also behaves in a similar fashion. 

Before proceeding further let us offer some comments on the operators that we considered. 
We have studied three different types of non-local operators as probes of thermalization. It is 
not a priori clear which of these operators provides the correct time-scale for thermalization. It 
can be explicitly checked (like in the case of vanishing chemical potential) that it is actually the 
entanglement entropy which thermalizes the latest and thereby sets the equilibration time-scale. 



This is demonstrated in fig. 25 
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Figure 25: For fixed (47rT)£ = 0.7 we have shown how Tent depends on X(4) for various probes. 
The red curve corresponds to the 2-pt function, the blue curve corresponds to the Wilson loop 
and the black curve corresponds to the entanglement entropy. Here T^^t corresponds to that of 
the two-point function. 



In the presence of a chemical potential the physics becomes richer. First, the relevant quantities 
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to consider here are the following dimensionless combinations: (Trcrit), (Ti) and (fJ^/T) — where 
T is the temperature of the thermal background, i is the length of the non-local operator, /i is the 
chemical potential and Tcru is the thermalization time. Thus, in general, (Trent) is a function of 
the two independent variables: {T£) and (fJ^/T) and our goal is to explore this behavior. 

Based on general physics intuition, once the thermalization has set in, we can identify the 
following different regimes 



Ti fixed, small : 



/i/T > 1 and fi/T < 1 
quantum classical 



(5.100) 



and 



/i/T fixed, small : and T£>^ . 

quantum classical 



(5.101) 



In the previous subsections, by analyzing the various probes of thermalization, we have seen 
qualitatively different features in precisely these "classical" and "quantum" regimes as defined 
above. The ubiquitous property that reveals itself through this exercise is that thermalization 
becomes faster for small values of /i/T, but for large values of /i/T the thermalization time in- 
creases without any upper bound. This non-monotonic behavior in the two regimes are smoothly 
connected through a minima as depicted in e.g. figs. [l7|(b), [22|c). This non-monotonicity is en- 
hanced in higher dimensions]^ We should also note that the identification of the "classical" or 
the "quantum" regime is only meaningful when the system has completely thermalized. Although 
it seems that the thermalization process is sensitive about this equilibrium feature, it is not clear 
to us whether such a "classical" or a "quantum" regime actually exist during the non-equilibrium 
period. 

Let us also briefly comment on the functional dependence of Trent on the dimensionless com- 
binations Ti and /i/T. If we take constant T£-slices, our generic observation suggests a linear 
relation 



Trcrit = A {d,Ti)^ for |»1 



(5.102) 



where A{d,Ti) is the slope of the curve which generally depends on the dimensionality of the 
problem and the value of Ti. On the other hand, if we take constant (/i/T)-slices we get the 



-"^•^ Although we have not considered the case of d = 2, which corresponds to an (1 + l)-dim dual CFT, we have 
explicitly checked that the non-monotonic behavior is further suppressed for d — 2. For small values of n/T, the 
variation in Trcrit is negligible and within the numerical accuracy it is difficult to conclude whether thermalization 
time increases or decreases. For larger values of /i/T, the dependence is similar to what we see in higher dimensional 
examples and Trcrit increases linearly with /i/T. 
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following polynomial relatioE 



14 



Tn,it = B {d, fi/T) {Tif + C {d, fi/T) (Ti) , (5.103) 

where B and C are two constants which depend on the dimensionality of the problem and the 
value of the chemical potential. Furthermore, it can also be checked explicitly that B{d,fi/T) < 
for fi/T <^ 1 and B{d,fi/T) > for /i/T 3> 1; whereas C{d,^/T) is always positive. For small 
values of the length of the operator, thermalization time behaves linearly with the length of the 
operator and the slope is given by C. For large values of the length, Tcrit is either sub-linear or 
super-linear depending on the sign of the other constant B. 

We have shown that the presence of a chemical potential makes the physics of holographic 
thermalization richer. Some important qualitative differences appear, namely, the non-monotonic 
behavior of the thermalization time as a function of /i/T . Some other charcteristics like top-down 
thermalization (UV modes thermalize first) and the fact that it is the entanglement entropy that 
sets the thermalization scale are common to the /i = case. 

This perhaps is a tip of the iceberg and there is a richer story waiting to be unfolded. Let us 
point out some open problems and future directions which would improve the understanding of 
holographic thermalization of theories with non-zero chemical potential: 

• Fortune favors the brave, or does it: The time-dependent background we have worked on is 
more "phenomenologically motivated" than one obtained from a "first principle" calculation 
even from the point of view of classical gravity. One might therefore wonder about the 
validity or the applicability of the physics we observe within this framework. In this work we 
have implicitly assumed the existence of at least some approximation in which the physics we 
see holds true. In [3] the authors showed that a dilaton source of small amplitude and finite 
duration at the boundary produces a wave that propagates in the bulk and collapses to form 
an AdS-Schwarzschild geometry. The spacetime for this collapse process was constructed 
as an expansion in the amplitude of the dilaton source and it was shown that it takes the 
Vaidya form at the leading order. A fascinating question on its own right is the study of 
the gravitational collapse process that may give rise to the AdS-RN- Vaidya metric. Perhaps 
a natural point to start from is to consider a generalization of [3] by considering a charged 
dilaton field. Such an exercise will either establish the AdS-RN- Vaidya background on a 
firm footing or we will learn some other interesting lesson. 

• The fact that thermalization time increases with increasing chemical potential (at least 
for large enough values of chemical potential as we have observed) can be reconciled with 

^''^This can be verified by explicitly fitting the data. 
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weak-coupling field theory intuition. For a bosonic system, increasing chemical potential 
"enhances" the Bose-Einstein condensation and therefore hinders the process of thermal- 
ization which will populate excited states. For a system with fermions, increasing chemical 
potential increases the available states by increasing the Fermi energy. Thus it is intuitive to 
conclude that a system with bosonic or fermionic degrees of freedom will thermalize slower if 



the chemical potential is increased^ However, to the best of our knowledge, we are unaware 
of a field theory (even toy-model) calculation that demonstrates this physics explicitly. On 
the other hand, for small chemical potential the faster thermalization that we observed does 
not fit our naive intuition. Thus a field theory computation at weak coupling will shed much 
light on the underlying physics and we may be able to isolate the features governed by strong 
coupling and the physics that is governed by the presence of a chemical potential. 

• It would be interesting to study the spread and evolution of correlations in an out-of- 
equilibrium system which eventually thermalizes. In [31] the authors pointed out that 
mutual and tripartite information are adequate probes to study this problem and calcu- 
lated these quantities holographically for a three dimensional bulk theory which is dual to 
an (1 + l)-dim CFT. Field theory computations of such observables are so far available only 
for the (1 -|- l)-dim. It will be an interesting exercise to generalize these computations in 
higher dimensions and analyze the corresponding physics. Furthermore, we have shown that 
a non-zero chemical potential introduces non-trivial new behavior in the thermalization of 
entanglement entropy. Thus, we expect that the mutual and tripartite information results 
will also be substantially modified when 7^ 0. 

• More phenomenological quantities can also be studied. For example the stopping distance of 
a massless particle — related in the weak coupling to the jet quenching parameter for QCD- 
like theories — was calculated in [32] for an AdSs-Schwarzchild black hole. An interesting 
question is how the stopping distance is modified in a time dependent background undergoing 
thermalization with and without chemical potential [36]. 

We hope to address some of these problems in the near future. 



6 Acknowledgements 

We would like to thank Willy Fischler and Berndt Miiller for conversations and encouragement 
about this work. E.G. acknowledges support of CONACyT grant CB-2008-01- 104649 and CONA- 
CyTs High Energy Physics Network. AK would like to thank the hospitality of the KITP, Santa 

^^We would like to thank Berndt Miiller for suggesting this possibility. 



39 



Barbara during the workshop "Novel Numerical Methods for Strongly Coupled Quantum Field 
Theory and Quantum Gravity" and financial support in part by the National Science Foundation 
under Grant No. PHYll-25915. AK is supported by the Simons postdoctoral fellowship awarded 
by the Simons Foundation. This material is based upon work supported by the National Science 
Foundation under Grant PHY-0969020. 



Appendix A. Geodesies in the background 

For the sake of completeness, let us first write down the AdS-RN-Vaidya background in general 
{d + l)-bulk dimensions (with d > 2) 

ds^ = — {-f{z, v)dv^ - 2dvdz + dx^) , = q{v)z'^-^ , (A.l) 

which gives 

2«:T- = (^-^z'-'^ - ^z^'-'g^) , ^JL. = {d- 2)L'-'^S^^ . (A.3) 



dv dv J dv 

Analyzing the general behavior of the geodesies or the minimal surfaces of various dimensions 
in the AdS-RN-Vaidya background is an interesting problem]^ For the rectangular cases we have 
considered in the main text, the equations of motion of such minimal surfaces take the following 
form 

zv" + {2p)z'v' -p + v'^ (^pf -^z^pj = 0, (A.4) 

where 

p = 1 =^ geodesic , (A. 5) 

p = 2 =^ Wilson loop , (A. 6) 

p = 3 =^ Entanglement entropy . (A. 7) 

On the other hand, the action functional for the minimal surface for the circular region is 

S= f dp^^/l- fv'^ -2z'v' . (A.8) 



""^^For a recent exhaustive study of such n-dimensional extremal surfaces in general static asymptotically AdS 
background of various dimensions, see e.g. [37]. 
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The equations of motion obtained for the circular region are involved, hence we do not present 
their explicit form here. We will briefly comment on a few general properties of the geodesies or 



the minimal surfaces obtained as solutions of (A.4). 



Before doing so, let us introduce the notion of an apparent horizon which will be relevant for 
the time-dependent backgrounds we have considered in the main text. A trapped surface T is 
defined as a co-dimension two spacelike submanifold with the property that the expansion of both 
"ingoing" and "outgoing" future directed null geodesies orthogonal to T is everywhere negative. 
The boundary of the trapped surfaces associated with a given foliation is defined as the apparent 
horizon. 



In what follows, we will closely follow [38]. For the background in (A.l) the vectors tangent to 
the ingoing and outgoing null geodesies are given by 

2 2 

/_ = -a. , l^ = -^d. + ^fd, (A.9) 

such that 

/_-/_=0, /+-/+ = 0, /_■/+ = -!. (A.IO) 

Now the volume element of the co-dimension two spacelike surface (orthogonal to the above null 
geodesies) is given by 

E^0)^". ,A.U, 

The expansions are defined to be 

^± = £±logS = /^9^(logS) , (A. 12) 

where C± denotes the Lie derivatives along the null vectors l±. The apparent horizon is then 
obtained by solving the equation B = 0, where 6 = 6'+6'_ is the invariant quantity. In this case 
we find 

= ^^^/ = O (A.13) 

gives the location of the apparent horizon which we can find numerically. 

Now we are ready to present a few representative figures demonstrating how a geodesic or a 
minimal area surface looks like in our time-dependent backgrounds. Such representative profiles are 



shown in fig. |26j It is clear that depending on the boundary separation length i, the corresponding 
minimal surface may or may not penetrate the apparent horizon. If the corresponding geodesic 
or the minimal surface penetrates the apparent horizon, the derivative of the curve becomes 
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discontinuous. The curve consists of two parts, one inside the apparent horizon and the other 
outside of it. These two sections of the curve are connected according to an analogue of the 
Snell's law of refraction [5]. This discontinuity is visible from the plots in the {z{x) — x} and the 
{z{x) — v{x)} planes. 



z(x) 
3.0 

2.5 

2.0 

1.5 

1.0 

0.5 

0.0 



" 



0.2 0.4 



0.6 0.8 

(a) 

z(x)-i 
3.0 

2.5 

2.0 

1.5 

1.0 



1.0 



1.2 



0.5 



v(x) 
1.5 

1.0 

0.5 

0.0 



I 
I 
f 
I 



0.2 



0.4 



I 
I 



0.6 

(b) 



0.8 



1.0 



1.2 



I 
I 
I 
I 



I 
I 

/ 
/ 



-0.5 0.0 



0.5 1.0 

(c) 



v(x) 



Figure 26: We have presented here the profiles of rectangular Wilson loop for d = 3. The black 



dashed line denotes the location of the apparent horizon obtained from solving equation (A. 13). 



We have presented here the profiles for only x > branch since the configuration has an exact 
symmetry under x —x. The blue, green and red curves correspond to (47rT)£ = 5.6, 8, 10.4 
respectively. 



An interesting feature that we have observed in e.g. fig. 21 and was also pointed out in [551 15] 
is the swallow-tail behavior of the corresponding evolution function. This is a generic feature that 
is observed from computing the various probes of thermalization by using solutions obtained from 
solving the equations in (A. 4) or the equations obtained from minimizing the action functional in 
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equation (A. 8) for large enough boundary separation. This swallow-tail behavior results from a 



multi-valuedness in as a function of the boundary time t as demonstrated in Fig 27 As has 
been pointed out in [5], this feature does not imply anything unphysical; we just need to be careful 
and follow the steepest descent procedure. 
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Figure 27: The multi-valuedness of with as a function of the boundary time t for d = 4 
rectangular entanglement entropy computation. The boundary separation length is {4irT)i « 2. 
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Appendix B. Charged AdS black holes 



Here we will briefly comment on some properties of charged AdS black holes in four and five bulk 
dimensions and their embedding in 10 or 11- dimensional supergravity. Let us begin with the case 
d = 4, i.e. in five bulk dimensions. 

The relevant theory here is the reduction of type IIB supergravity [391 HO] truncated to the 
M = 2 sector with an U{lY symmetry. The corresponding five dimensional theory [JT] includes 
three U{1) gauge fields, denoted by Ai (/ = 1, 2, 3), and two scalar fields, denoted by X/ subject 
the constraint X1X2X3 = 1. This theory has a family of black hole solutions parametrized by 
three charges (corresponding to the three f/(l)s) and a non-extremality parameter usually denoted 
by n [221 [20] The general form of such black hole solutions take the following form 

ds^ = {H^H^H.r'/' fdt' + {H,H2H,Y^' (^r'dQl, + ^^ , (B.l) 

2 

Hi = 1 + ^, f = k-^+'^ (HiH^Hs) , Aj = [HJ^ - l) dt , (B.2) 

where k takes values or 1 depending on whether we work in the Poincare patch or the global 
patch: the line element denoted by dQf q in this case corresponds to that of M.^. The three charges 
qi, when uplifted to 10-dim type IIB supergravity, correspond to angular momenta along three 
Cartan directions on the five-sphere. We have considered only the Poincare patch solution in the 



main text. Furthermore, the AdS-RN solution that we have considered in e.g. equation (2.4) 
corresponds to the case when qi = q2 = Qs = q- In this case, the scalars become trivial and the 
background metric takes a more familiar form 

- ^) , (B.4) 

where we have defined = r"^ + q and ph is the location of the event-horizon. We can obtain the 



background written in equation (2.4) from the above expression by identifying the following 



M^A, ,B.5) 

It was noted in |PU], in order to have a horizon the non-extremality parameter p must satisfy 
a lower bound which can be analytically obtained. This means for a given /i, there is an upper 
bound for the charge beyond which there is no horizon. This fact is revealed in the main text in 



"'^''This non-extremality parameter /i is not to be confused with the chemical potential that we have defined in 
the main text. 
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e.g. fig. [T| If we violate tliis bound, tlien a naked singularity appears |20]. As shown in ^42j, this 
naked singularity has a natural meaning within string theory as an ensemble of giant gravitons 
which are distributed over the compact S^. For more details, we will refer the reader to [H]. Note, 
however, that even within this bound the relevant dimensionless quantity X{d), defined in equation 



(3.22), has a range given by X{d) ^ [0,oo]. Thus the existence of an upper bound in the charge 
Q does not imply any corresponding bound for the chemical potential measured in units of an 
appropriate power of the temperature. 

It is transparent that if the charge violates this upper bound, then the corresponding physics 
does not describe a thermal state in the dual field theory. Our primary focus in this article was 
to study the thermalization process by analyzing various probes of thermalization. It is intriguing 
to ask what will happen to these probes if we carefully fine-tune the mass parameter M and the 
charge parameter Q such that the background cannot form a horizon. In this case, the probes 
should never thermalize. Following the procedure explained in the main text, we can obtain the 
behavior of probes of thermalization in such a case and one such representative plot is shown in 



fig. 28 
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Figure 28: The behavior of Ahermai (measured in units of the AdS-radius) as a function of 
boundary time for d = 4. We have set the boundary separation length ^ = 2 in units of the black 
hole mass. The blue curve corresponds to Q = 1 and the red curve corresponds to Q = 2 (in 
units of the black hole mass). Note that with these values of the parameters, the corresponding 
AdS-RN background does not have an event-horizon; instead it develops a naked singularity. 



It is clear from fig. 28 that initially Ahcrm begins to follow a similar bahviour as in the case 
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when they eventually thermalize; however, after a certain amount of time has elapsed, the curve 
turns back and never really achieves thermal equilibrium. For larger values of Q (with M fixed), 
this "turn over" behavior is further enhanced. Note that in the absence of a horizon it is not 
entirely clear to us what the corresponding non-local observables correspond to; hence we just 
present one representative plot and refrain from any further claims. 

Similarly, it is also possible to consider the S* ''-reduction of 11-dimensional supergravity which 
gives rise to 50(8) gauged A/" = 8 supergravity in four dimensions. As in the five dimensional case, 
it is possible to take an A/" = 2 consistent truncation which gives rise to a bosonic sector consisting 
of a graviton, four Abelian gauge potentials, three axions and three dilatons. This theory also has 
charged black hole solutions parametrized by four charge parameters and is given by 

ds^ = -{H,H,HsH,)-'^'fdt' + {H,H2HsH,f'(^r'dnl, + ^^ , (B.6) 

2 

Hi = 1 + ^ , f = k-^ + ^ (H.H^H.H,) , Aj = {Hj^ ~ l) dt , (B.7) 
I = (B.8) 

This class of charged black holes can be embedded in M-theory in terms of rotating M2-branes. 
It is again straightforward to check that the usual AdS-RN background is recovered with four 
charges set equal. 
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